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Abstract 

I present a calculation of the interfacial width within the capillary wave 
(Gaussian) approximation. The calculation is done on rectangular 
lattices of size Li times L2, with periodic boundary conditions. 



1 Introduction 



The two-dimensional interface separating the phases of a 3D binary system 
can undergo a roughening transition Below the roughening temperature, 
the interface has a finite width in the infinite volume limit. The rough phase, 
however, is characterized by strong fluctuations of the interface position vari- 
ables, and the interface width diverges when the area of the surface goes to 
infinity. Properties of the roughening transition have been investigated in a 
variety of studies, see e.g. and references therein. It is generally believed, 
and has been confirmed in a number of numerical studies 0, that the in- 
frared (large distance) properties of a rough interface can be described by 
massless Gaussian modes. This concept was introduced by Buff et al. |^ and 
is usually referred to as the capillary wave model (CWM). The CWM allows 
to make predictions, e.g., about the finite size effects of various interface 
properties ||, §. Quite some attention has been devoted to asymmetric in- 
terfaces, where the two extensions Li and L2 of a rectangular interface do not 
necessarily coincide ||, H, |10|. In [H, ^ an extended (non- Gaussian) 



CWM was employed with an action proportional to the interface surface. 
The interface shape dependence was investigated in a field theoretic setting 
in . 

In this article, I present a calculation of the interface width on asym- 
metric rectangular lattices, within the CWM. Related calculations with a 
continuum cutoff were done in a gauge theory setting in Interestingly, 
the lattice result obtained in the present paper points the way to a significant 



simplification of the continuum result presented in |Tl 



2 The Lattice Capillary Wave Model 

Let the interface (without overhangs) be described by a "height" function (p^, 
where x is a site of a two-dimensional grid with side lengths Li and L2. We 
choose units such that the lattice spacing is dimensionless and assume that it 
is one for both lattice directions. We employ periodic boundary conditions. 

Within the capillary wave approximation for the fluctuations of the height 
variables, expectation values are defined by 



1 



with the Hamiltonian (energy functional) 

Hoiv) = ^J2i'Px- vyf + ^ E 

{x,y) ^ 

The first term is a sum over all nearest neighbour bonds in the lattice. We 
define the interface width (or surface thickness) by 

where xq is an arbitrary site of the lattice. The interfacial width can also be 
written as 

= 2 ((^,„ - 0)2) , 

where 

X 

is the average interface position. In the following section we will show that 

^ + ^ ("'-''''(»)) + ^2M + o((i A)-) . (1) 

77 is the Dedekind function. We will compute explicit expressions for the 
constant K and the function Z{u) that depends on the asymmetry-parameter 

3 Calculation of W'^ 

By means of discrete Fourier transformation W"^ can easily shown to be 

2,3 1 



The lattice momenta are defined hy pi — rii, with = 0, . . . , Lj — 1, and 
pf = 4 sin^ Let us rewrite the sum excluding p — {pi,P2) = in the 
following way: 

/ L2-1 L1-1L2-1 

E /K' ^2) = E /(O' ^2) + E E /(^i' ^2) ■ 

n2=l ni=l n2=0 



Employing the identity 
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E^—x\n\—iujn 




sinh (x/2) + sin (u;/2) sinhx 
we can perform a first summation: 

sin^ (Hii ) + sin^ (h^) sinh X^, [ 1 - e'^^x^, 
The quantity is defined by 

smh = sm 



Performing a suitable limit, one can use eq. @) to show that 



L2-1 I I 



E — ^ = -{Ll~l). 



n2=i sm 

Putting things together, we find 

'2^~V2 L1L2 ^ 2Li „f^, sinhX„, 1^ ^ 1 - e'^^x^, 
We shall first study the sum 

1 ^^"^ 1 
S = — — V 

2Li„f^, sinhX„, ■ 

Note that 



sinhX„, = 2 sinh l'^) cosh ( ^] = 2 sin f ^) Wl + sin^ 



Exploiting that sin (vmi/Li) = sin {tt{Li — rii)/ Li) and assuming that 
2Mi is even, we obtain 

Mi-l 

5 = 2 5: ^?(X„J+^(MO, 



with g{xni) = (2Li sinh(X„J) ^. S can then be represented as S = Si + S2, 
where 



A/1-1 / 

ni=l 



1 



2Li sin (f^) Wl + sin^ (nii) 27^1 



1 



AV2L1 



and 



The sum 5*1 can be evaluated with the help of the Euler-Mc-Laurin summa- 
tion formula. One finds 



27T \2 

5*2 can also be summed 



S, = — {C + i^{M,)) , 

where C = 0.5772... denotes Euler's constant, and ip is the Digamma func- 
tion. We expand 
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S thus combines to 



-ln2 + lnLi + C 



1 



1 1 



Li 3Li 



5 = — C + -ln2-ln7r + lnLi + ^^ZTo + O [Ll 



2tx 



Our next task is to study 



TT 



1 



72L2 



-4 



Li sinh(X„J 1 - e-^2X„, • 
As for eq. (||) , we use the symmetry of the "integrand" to write 



Afi-l 



Ml J 



ni=l 



(4) 



with 



Li sinh(X„J 1 - e'^'^^^i ' 

Recall that Mi = Li/2, and that we have assumed that Li is even. We 
are interested in the limit Li — > oo, with the ratio u = L2/L1 kept fixed. 
One first observes that the contribution h{XMi) vanishes exponentially fast 
in this limit, namely like exp(— 2 arsinh(l) -uLi). For large Li, we expand 

1 o'^i 

M^»i) = - ^^(l _gn,) il + A{ni,U,u)) , 

where 

q = exp{—2nu) . 

We obtain 

Extending the sum in eq. to the range 1 to infinity introduces only errors 
that are exponentially small in Li. We thus find 

vr„^,n(l-g-) + L? I 3 (I - Q^)' 3 „^ 1 - a" i ^ ^ ^ 

Let us define 



[~ 


00 

n=l 


F{u) 


= E: 


G{u) 


= E: 


H{u 


)=E: 



n{l-q") ' 



Noting that InLi = ^ ln(LiL2) — | Inw, and putting everything together, we 
obtain 



\nL1L2 1 

-7777 = K + 



- In (wV4 ^(,^)) + -^Z{u) + O . (5) 



2(3 47r TT ^ " LiL 

The constant i^T is given by 

ir = — fc + iln2-ln7r 
27r V 2 



Furthermore, we have used that F{u) = Gi{q), where 



; nP (1 - x"") 



is Lambert's series [0|. For p = 1, it is related to Dedekind's r^-function, 

G, (e-^-")=-^-lnr^M, 

with 

oo 
n=l 

The symmetry of the interface width under exchange of Li and L2 should 
be reflected in an invariance under the transformation u 1/u. It is well 
known that u^/^rfiiu) is invariant under this transformation. 
The part proportional to {LiL2)~^ is given by 

Ziu) = —u- — + —u" Giu) - - uHiu) . 
^ ' 11 12 3 ^ ^ 3 ^ ' 

We would like to demonstrate that this expression is invariant under u ^Xju. 
To this end we first observe that 

i7(n) = ^(l-E2M) , 

with E2 being the first Eisenstein series. It obeys the functional relation 

6 

E2{iu) = -u-'^E2{iu'^) + -u-\ (6) 

TT 

The crucial step is now to recognize that G can be written as 
We then have 

Z{u) = ^ (uE2{tu) + u'^E2{zu)^ - ^ . (7) 

Differentiating eq. (P) with respect to u yields the behaviour of -^E2{iu) un- 
der u ^ 1/u. Using this and eq. (|^), it is easy to demonstrate the invariance 
of the right hand side of eq. (^. 
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We remark that it would be easy (though technical) to extend eq. to 
higher orders in (LiL2)~^. 

A comparison with eq. (A. 4) of ref. [jlT| leads to an interesting observa- 
tion. For the interface thickness on the continuum torus, regularized by a 
point splitting procedure, the authors obtain (I adapted their notation to the 
present setting) 



3 7TU 

2 12 



2'K(tW'^ = In (LiVI + uVSe) + ^ arctan u + | arctan ^ 

- ^k=i ^'^^S:^ ((1 + u>'t B', sin m + 2) arctan t.) . (8) 

Here, 



is the fc'th Eisenstein series, and -B^ are Bernoulli numbers, defined through 

j;^-i-i-gMf(|-.-. (9) 

Comparing eq. (||) with the infinite area limit of eq. (|^), and identifying 
a = 1/(2/9), one is led to conjecture the identity 

^\n{l + u'^) +g{u) = -2\n7]{iu) -C , (10) 

with 

g{u) = 2^ arctan M + | arctan ;^ — | — ff 

- ^k=i ^^SS&ir ((1 + «')^')' B', sin ((2fc + 2) arctan^) 

(11) 

Eq. (|TUp turns out to be true, with C = —ln{TT), thus leading to an enormous 
simplification of the result of |Tl| . 1 will only give a sketch of the proof. One 

uses that 

sin((2A; + 2) arctan(M)) = (1 + M^)"^'"^Im (l + iuf'^'^ 
One then meets the following two expressions to be evaluated: 

-| oo p/ 
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and 



V 



-Im -—- — — (ravrfl + iu)Y 



7r% ^,n3(l-g")^tt(2A; + 2)! 

The latter expression can be evaluated by first doing the fc-sum (yielding a 
cos minus the two leading terms) and then extracting the imaginary part. 
One then again recognizes Lamberts series with p = 1, and ends with 

i\u 

V = 2 ln?7(m) . 

6 

For the study of Y, formula (A. 6) of [|11| is very helpful: 

00 0/ oo -mz _ 1 ^2 2 q 3 

2k{2k + 2)\' m3 ^ 2^""'^ 6 12- 

With z = 7t{—u + i), and noting that for m > 

ln2; = Invr + - ln(l + u"^) — i arctan(l/M) + ire , 
one arrives after some algebra at 

1 1 1 /i 2\ 1 I U 1 3 TTU TC 

r= — Invr ml + w H arctan arctan — I 1 . 

2 ^ ^ 2u u 2 u 2 A 4u 

Combining things and noting that (for u > 0) arctan u + arctan ^ = 7r/2, one 
arrives at eq. (p!OD. 



4 Concluding Remarks 

The calculation presented in this paper is interesting at least for five reasons. 
First, it constitutes an illustrative example of analytical techniques that allow 
the exact evaluation of certain lattice sums. Second, the invariance under 
the exchange of the two lattice directions leads necessarily to the occurrence 
of modular forms and interesting relations between them. Third, the lattice 



result helped to simplify very much the continuum computation of ref. |11 



Fourth, interesting physical applications are possible, e.g., in the context of 



analytical and Monte Carlo calculations (cf. again JTTl). Finally, the present 
approach gives the the 1/area correction to the asymptotic behaviour of the 
interface thickness. This could be useful as a starting point to test the CWM 
"beyond the gaussian approximation" P, 110] also for the interfacial width. 
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